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Abstract 


Navar, J., € Lizárraga-Mendiola, L. (January-February, 2016). 
Temporal Patterns in River Flows in Basins in Northern 
Mexico. Water Technology and Sciences (in Spanish), 7(1), 35- 
44. 


Tn the dry basins in northern Mexico, concern about rivers is 
focused on seasonal and long-term variations in river flows 
that may be attributable to climate change, since it affects the 
water supply for domestic, industrial, public and agricultural 
uses. The objective of this research was to separate daily 
flow components and understand the patterns in 172 time- 
series corresponding to five rivers in northern Mexico. The 
headwaters of the San Pedro, Sinaloa, Nazas-Aguanaval, San 
Juan and San Fernando-Soto La Marina rivers are located in 
the Sierra Madre Occidental and Oriental, and flow through 
the plains into the Pacific Ocean, the Gulf of Mexico and 
the Central Plateau in northern Mexico, providing over 
90% of the water supply for all uses. The Rora computer 
program isolated the base and direct flows from the total. 
The Mann-Kendall test and the Seno S method evaluated 
the trends and magnitudes of the changes in annual and 
monthly flow, monthly standard deviation, and daily 
maximums and minimums for total, direct and base flows. 
The results found that the base and direct flows were in 
equilibrium for most of the rivers. Meanwhile, each one of 
the variables analyzed had statistically significant trends at 
40% of the hydrometric stations. Over 26% of the stations 
showed a consistent reduction in flow over time. Although 
additional investigations are needed to quantitatively 
isolate the potential sources of the variations to explain the 
temporal trends, the causes can preliminarily be attributed 
to climate variability in combination with natural resource 
management practices. The effect of potential climate change 
cannot be excluded from these patterns. 
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Resumen 


Navar, J., € Lizárraga-Mendiola, L. (enero-febrero, 2016). Patrones 
temporales de flujo del río en las cuencas del norte de México. 
Tecnología y Ciencias del Agua, 7(1), 35-44. 


En las cuencas áridas del norte de México, las preocupaciones por los 
ríos se centran en las variaciones estacionales y de largo tiempo en 
la descarga de los ríos que pudiera ser atribuida al cambio climático 
porque controla el suministro de agua para los usos domésticos, 
industriales, públicos y agrícolas. El objetivo de esta investigación 
fue separar los componentes del caudal diario y entender los 
patrones de 172 series de tiempo de cinco ríos del norte de México. 
Los ríos San Pedro, Sinaloa, Nazas-Aguanaval, San Juan y San 
Fernando-Soto La Marina nacen en las montañas de las Sierras 
Madre Oriental y Occidental, y fluyen a través de las planicies del 
océano Pacífico, el Golfo de México y la Mesa Central del norte de 
México, proveyendo con más de 90% del suministro de agua para 
todos los usos. El programa de cómputo Rora separó los flujos 
base y directo del total. La prueba de Mann-Kendall y el método 
del Seno S evaluaron las tendencias y magnitud del cambio de la 
descarga anual, mensual, desviación estándar mensual, y mínimo 
y máximo diario para los flujos total, directo y base. Los resultados 
mostraron que los flujos base y directo se encuentran en equilibrio 
para la mayoría de los ríos, pero más de 40% de las estaciones 
hidrométricas tuvieron tendencias estadísticamente significativas en 
cada una de las variables analizadas. Más de 26% de la estaciones 
mostraron un reducción consistente del caudal en tiempo. Aunque se 
necesita de investigación adicional para aislar de modo cuantitativo 
las potenciales fuentes de variación que explican las tendencias 
temporales, éstas pueden ser preliminarmente atribuidas a la 
variabilidad climática, en combinación con las prácticas de recursos 
naturales, pero el efecto del potencial cambio climático no puede ser 
excluido de estos patrones. 


Palabras clave: tendencia negativa, magnificación de la sequía, 
caudal base, directo y total. 
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Introduction 


Conventional management practices of water 
resources cannot meet future water demands 
for development in arid, semi-arid, and sub- 
tropical countries (UNCED, 1992; Postel, 2000; 
Comprehensive Assessment Secretariat, 2006). 
Mexico's northern watersheds span the arid 
and semi-arid belts of North America of the 
boreal Hemisphere; where river flow is featured 
by erratic and variable temporal and spatial 
discharge (Ward, 1998; Návar-Cháidez, 2012). 
In addition, agriculture, population, industry, 
trade and services, among others rely upon 
river flow to meet conventional water demands 
(Conagua, 2012). In Mexico's northern drylands, 
these economic sectors have been steadily grow- 
ing over the past 60 years due to (a) population 
growth and (b) economic development by in- 
creasing commercial treaties with other nations. 

Climate variability coupled with steady 
increasing water diversions modifies river flow 
and takes a toll on wetland ecosystems (Postel, 
2000; Návar-Cháidez éz Lizárraga-Mendiola, 
2013). The perturbation of river flow then re- 
quires new strategies to cope with conventional 
and environmental water demands, creating 
a positive feedback loop between increasing 
river perturbation and need for adaptation. For 
several watersheds, during acute dry episodes, 
river flow does not meet conventional demands 
and therefore society currently rely upon deep 
drilling and the inter-basin transfer of water 
(Conagua, 2012). In parallel to inter-basin water 
transfers, several sustainable practices of hydro- 
logic resources are occurring (Conagua, 2005). 
However, unsustainable practices remain on 
these watersheds, such as the high volume of ir- 
rigation vs low productivity of several irrigation 
districts; under valuing water resources, and 
highly regulating flow of most rivers (Pereyra, 
Oweis, €: Zairi, 2002; Comprehensive Assess- 
ment Secretariat, 2006). In some cases, water use 
efficiency has increased, resulting in doubling 
productivity with the same amount of water. In 
addition, water withdrawals from domestic use 
(sometimes referred to as grey water) are being 


+» ISSN 0187-8336 


treated to further irrigate agricultural fields, 
along with other sustainable practices (Postel, 
2000). 

Climate variability is responsible for er- 
ratic and infrequent high rainfall episodes and 
river flow anomalies that usually start or end 
long drought episodes, a common feature of 
Mexico's northern as well as of most semi-arid 
drylands. Hurricanes and cyclones entering 
northeastern Mexico have a frequency period of 
close to 15 years (Návar-Chaidez, 2012). These 
climatic phenomena restore aquifer and reser- 
voir storage, as well as river flow components. 
In dry spells, storage is frequently depleted. 
Drought episodes of several temporal scales; 
inter-annual, 3-7 years, 9-12 years, and 60-80 
years have been well documented in the region 
and they have been correlated with large-scale 
climate phenomena such as El Niño /Southern 
Oscillation, the Pacific Decadal Oscillation, the 
Atlantic Multidecadal Oscillation, or by a com- 
bination of any of these synoptic-scale climate 
events (Cavazos éz Hastenrath, 1990; Comrie éz 
Glenn, 1998). Climate change may also in turn 
contribute to magnify river flow variability of 
most arid and semi arid watersheds. Accord- 
ing to Mulholland et al. (1997), whose projec- 
tions are based on Global Circulation Models, 
Northern Mexico would likely receive between 
5 to 20% less rainfall and between 5 and 25% 
less streamflow should global warming trends 
continue. 

The decadal drought spells of the 1950's, 
the 1990's, and the inter-annual droughts of the 
2010's have prompted renewed efforts to quan- 
tify water storage and the length of temporal 
drying patterns. Hence, this report addresses 
whether river flow components of 172 gaug- 
ing stations from five rivers show temporal 
monotonic patterns in the hydrologic regions 
of Sinaloa (RH10), Presidio-San Pedro (RH 11), 
Nazas-Aguanaval (RH 36), San Juan (RH 37), 
and San Fernando-Soto La Marina (RH 25). 
The working hypothesis was that there is no 
consistent change in any of the streamflow pa- 
rameters or components over time and therefore 
no change in the magnitude of discharge. 
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Materials and Methods 
The Study Area 


The study area covers the basins of Sinaloa, 
Presidio-San Pedro, Nazas-Aguanaval, San 
Fernando-Soto La Marina and the Rio San Juan. 
These watersheds span the Mexican northern 
States of Sinaloa, Nayarit, Durango, Chihuahua, 
Zacatecas, Coahuila, Nuevo Leon, and Tamau- 
lipas. 

Present and future demographic information 
on each watershed is reported in Table 1. 

A high diversity of climates features North- 
ern Mexico (Conagua, 2005). Cold climate is 
typical of the highest peaks and the uplands of 
the main mountain ranges. Temperate climates 
characterize the central Valleys and Mesas of the 
Cordilleras. Subtropical and tropical climates 
facet the low ranges of the plains of the Pacific 
Ocean and Gulf of Mexico. Dry climates are 
typical in the Central Chihuahuan Desert and 
the southern Great Plains of North America. 

Pine, oak, and a mixture of pine-oak forests 
cover extensive upland areas of the eastern 
and western Sierras Madre mountain ranges. 
Montane forests, dominated by thorny shrubs 
and low, broadleaved trees, distribute widely 
on the landscapes of the piedmonts of the main 
mountain ranges. Tropical and subtropical 
flora distribute in the plains and low lands of 
the piedmonts facing the oceans. Chaparral, 
dry, xerophytic plant cover is typical at the 
central Chihuahuan Desert and the southern 
Great Plains of North America. Tamaulipan 


thornscrub, acacia and mesquite forests extend 
into the lowlands of the Pacific Ocean, northern 
Gulf of Mexico and the Great Plains of North 
America (Rzedowskii, 1980). 

Shallow Litosols are typical of the main 
mountain ranges. Deep Regosols, Xerosols, 
and Yermosols characterize the landscape of the 
central Chihuahuan Desert. Deep Vertisols are 
distributed in the plains in the lowlands of the 
Pacific Ocean and the northern Gulf of México. 


Methodology 


Atotal of 172 hydrometric stations placed across 
five major rivers were selected for breaking 
streamflow components and testing the time se- 
ries for temporal monotonic patterns (Figure 1). 
The number of hydrometric stations analyzed 
per basin is depicted in Table 2. 


Hydrograph Separation Analysis 


The separation of multiple river hydrographs 
into base and direct flow is a difficult task. To 
simplify this assignment for long term river flow 
time series, total daily discharge was split into 
base and direct or quick flow using the RORA 
hydrograph separation techniques (Rutledge, 
1998). Before daily discharge data was fed into 
the RORA computer program, it was checked 
for: a) the presence of long periods of recession 
and b) the absence of intermediate missing data. 
Screen, a computer program, checks data for 
recession and continuity. The recession index 
is required in order to define the master curve. 


Table 1. Demographic growth of five Mexico's northern watersheds (Source: Conapo, 2000). 


Watershed Inhabitants 
2006 2015 2030 
1. Sinaloa (RH 10) 22791975 2 431 372 2612141 
2. Presidio-San Pedro (RH 11) 771 164 782 424 781 339 
3. Nazas-Aguanaval (RH 36) 1766 811 1868 697 1977 106 
4. San Juan (RH 37) 4 651 828 5 242 335 6 069 918 
5. San Fernando-Soto La Marina (RH 25) 682 033 727 322 781 397 
Total 10 151 811 11 052 150 12 221 901 
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Figure 1. Map depicting several hydrometric stations of five northern watersheds of Mexico. 


Table 2. The number of hydrometric stations analyzed for temporal patterns and magnitude of river discharge 


in five Mexico's northern watersheds. 


Watershed Hydrometric Stations a 
1. Sinaloa (RH 10) 58 1960-1999 
2. Presidio-San Pedro (RH 11) 15 1970-1998 
3. Nazas-Aguanaval (RH 36) 33 1940-1999 
4. San Juan (RH 37) 30 1955-1999 
5. San Fernando-Soto La Marina (RH 25) 36 1970-1999 


Data analyzed from Conagua (2005) records. 


Once the master curve is developed, PART is 
used as a module of RORA to separate base and 
direct flow. 

Monthly total discharge, standard deviation 
of daily discharge (within a given month), an- 
nual discharge, annual base flow, and annual 
direct flow or quick flow are parameters and 
components derived from daily discharge 
data. The nonparametric Mann-Kendall sta- 
tistic tested for temporal monotonic trends on 
the river flow time-series data. This test reveals 
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the presence or absence of temporal monotonic 
increasing or decreasing patterns. 

The non-parametric Sen's method quantified 
the magnitude of slope change in all five dis- 
charge variables. The slope coefficient describes 
a change of the variable per unit time (annual 
time scale). The Sen's method assumes a mono- 
tonic tendency in the time series. 

Data on land-use changes was compiled for 
several Mexico's northern ecosystems. In addi- 
tion changes in precipitation and temperature 
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parameters were also evaluated for several 
climatic stations. Correlations between river 
flow parameters and large-scale climate events 
(El Niño /Southern Oscillation, ENSO, The Pa- 
cific Decadal Oscillation, PDO, and The Atlantic 
Multidecadal Oscillation, AMO) previously 
reported aided to understand sources of river 
flow variability. This information, in addition 
to population growth, was linked to potential 
subtle river flow temporal patterns. 


Results 


Table 3 depicts the statistics of river flow com- 
ponents for all five studied watersheds. Base 
and direct flow are statistically similar for all 
watersheds. That is, on the average, 50% of river 
discharge is the low extended base flow rising, 
from shallow (mountain range uplands) and 
shallow bottom-land aquifers. The remaining 
50% of discharge is direct flow produced dur- 


ing and immediately after the rainy season that 
lasts from June to September. Hortonian surface 
and quick subsurface flow control direct flow. 
Watersheds draining into the Pacific Ocean 
(RH10 and RH11) and into the interior Central 
Plateau (RH 036) have a unimodal-monsoonal 
type of precipitation, with a major peak in July. 
On the other side, watersheds that drain into 
the Gulf of Mexico feature a bimodal type of 
rainfall distribution with two peaks, the first 
one during May-June and the last one during 
September-October. 

Table 4 depicts the number and percentage 
of gauging stations with statistically significant 
negative (-) and positive (+) temporal mono- 
tonic tendencies. 

In general, over 40% of each river flow pa- 
rameter analyzed showed statistical significant 
temporal monotonic trends. That is, the river 
flow time-series are not stationary in the first 
momentum. However, of the total, on the av- 


Table 3. Baseflow and quickflow statistics for several rivers of five Mexico's northern watersheds. 


Discharge Variables 
Watershed Baseflow Quickflow 
(+ confidence intervals) (%) (+ confidence intervals) (%) 
1. Sinaloa (RH 10) 0.52 (+ 0.08) 0.48 (+ 0.08) 
2. Presidio-San Pedro (RH 11) 0.52 (+ 0.08) 0.48 (+ 0.08) 
3. Nazas-Aguanaval (RH 36) 0.56 (+ 0.09) 0.44 (+ 0.09) 
4. San Juan (RH 37) 0.57 (+ 0.08) 0.43 (+ 0.09) 
5. San Fernando-Soto La Marina (RH 25) 0.50 (+ 0.06) 0.50 (+ 0.08) 


Table 4. Number of gauging stations with statistical significant temporal monotonic patterns of several river flow parameters. 


River Flow Parameters 

Watershed Annual Monthly Std. Dev. Base Flow Direct Flow 
Sinaloa 25(-4, +21) 26(-5, +21) 28(-17, +11) 28(-3, +25) 27(-16, +11) 
San Pedro 6(-4, +2) 5(3, +2) 6(-6, +0) 2(-0, +2) 6(-6, +0) 
Nazas 18(-15, +3) 18(-15, +3) 19(-16, +3) 11(77, +4) 18(-15, +3) 
San Juan 1311, +2) 14(-12, +2) 12(-9, +3) 12(8, +4) 12(-8, +4) 
San Fernando 12(-9, +3) 12(-9, +3) 8(-7, +1) 18(-14, +4) 8(7, +1) 
Total 74(-43, +31) 75(-44, +31) 73(-55, +18) 71(-32, +39) 71(-52, +19) 
Ratio (%) 43(25, 18) 44(26, 18) 42(32, 10) 41(19, 22) 41(30, 11) 


Std. Dev. = Standard deviation of monthly discharge. The last row was estimated from the ratio of the number of gauging stations with statistical 
significance to the total number of gauging stations analyzed. 
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erage, 26% had a negative and the remaining 
16% had a positive tendency. Total annual and 
total monthly river flow parameters are steadily 
declining at the San Pedro, Nazas, San Juan, and 
San Fernando watersheds but it is increasing at 
most gauging stations along the Sinaloa rivers. 

Most gauging stations with declining tempo- 
ral river flow patterns span the arid, semi-arid 
and subtropical watersheds RH11, RH36, RH37 
and RH25, characterized as having an important 
area in the interior Valleys of the Chihuahuan 
Desert and the southern Great Plains of North 
America. In these physiographic regions, mon- 
sonic erratic, infrequent rainfall and high evapo- 
transpiration losses control the hydrologic cycle. 
The human dominance on river flow is of least 
importance on these watersheds and therefore 
river flow is not as regulated asitis in e.g., the 
RH10. 

Most gauging stations with positive tem- 
poral monotonic trends in total annual and 
total monthly river flow are distributed along 
the Sinaloa (RH10) basin, a watershed with 
most hydrometric stations placed along the 
lowlands of the Pacific Ocean. Most gauging 
stations were placed below irrigation districts, 
reservoirs, and urban centers, where human 
control exerted on water resources (increased 
diversions from upland reservoirs as well as 
on deep groundwater resources and discharge 
of excess water in other places) may explain 
the positive temporal trend of river discharge. 
Therefore, increasing river flow patterns over 
time may be partially explained by water re- 
source management practices. This effect masks 
in the short run the potential subtle, changes in 
the local hydro — climate variables of the Pacific 
Ocean lowlands. 

Reduced direct flows and increased base 
flows over time could be explained by the 
regulated flow effect by man-made reservoirs, 
which reduce discharge variation of most gaug- 
ing stations below reservoirs when rivers are 
used as waterways to convey water to irrigation 
districts. Several river sections are experiencing 
the presence of smaller direct flows and larger 
base flows. Steady discharge of municipal sew- 


erages and irrigation excess into streams may 
partially explain these observations. 

Base flow is increasing in most gauging sta- 
tions along the Sinaloa and San Pedro water- 
sheds, but it is declining at most of the remain- 
ing gauging stations in the watersheds Nazas, 
San Juan, and San Fernando. Irrigation —excess, 
municipal discharge and increased diversions 
from deep aquifers of the interior watersheds 
partially explain this tendency. Subtle shifts in 
rainfall frequency, depth, and intensity may also 
be contributing to increasing aquifer recharge 
that may in turn control base flow discharge. 
In several watersheds, changes in the frequency 
and depth of rainfall coupled with increasing 
mean annual temperature may in turn contrib- 
ute to modify parameters of river discharge. 

Direct flows are also declining in most gaug- 
ing stations of all watersheds. Increased storage 
by man - made reservoirs above 30 000 Mm? 
in arid, semi-arid, and subtropical landscapes 
may partially explain the control exerted of 
precipitation-excess on river flow. The rate of 
change of discharge parameters is reported in 
Table 5. 

Most slope values statistically deviate from 
zero stressing the statistical significance on 
the declining of direct discharge over time. 
For total annual discharge, the San Fernando 
(-16.47 Mm? year”), Nazas (-4.74 Mm year”) and 
Sinaloa (-3.66 Mm? year*) watersheds have the 
largest reductions of total annual discharge over 
time. The remaining watersheds present total 
annual discharge of less than -1.00 Mm* year”. 
For monthly river flow, discharge reductions are 
more important at the watersheds San Pedro 
(2.05 Mn year”), San Juan (-2.10 Mm* year”), 
and San Fernando (-2.56 Mm* year”). 

Temporal patterns of peak flows had been 
notoriously modified (having on the average a 
reduction of 4.7 m* s”) in contrast to the small 
base flow reduction (average reduction of 1.4 
m? s*) for all watersheds. The Sinaloa and San 
Fernando basins present the largest human con- 
trol on river flow and therefore on peak flows. 
This, in turn, correlates well with the number of 
man-made reservoirs. Although base flow has 
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Table 5. The median slope (and confidence intervals) for five discharge parameters with negative temporal tendency of 172 


gauging stations placed in five Mexico's northern watersheds. 


River Flow Parameters (Mm year”) 

Annual Monthly Std. Dev. Base Flow Direct Flow 
Sinaloa -3.66(2.19) -0.82(1.06) -0.32(0.41) -1.91(0.51) -9.69(6.25) 
San Pedro -0.11(0.12) -2.05(0.94) -2.52(4.89) -1.02(0.77) 0(0) 
Nazas -4.74(3.98) -0.40(0.31) -0.50(0.28) -1.31(0.92) -1.28(0.80) 
San Juan -0.84(1.11) -2.10(0.92) -0.01(0.01) -1.58(0.61) -2.73(3.30) 
San Fernando -16.47(8.17) -2.56(2.03) -2.53(1.02) -0.12(0.10) -10.01(2.61) 
Average -5.16(3.11) -1.59(1.05) -1.18(1.32) -1.19(0.58) -4.74(2.59) 


Std. Dev. = Standard deviation of monthly discharge. 


been little modified in contrast to peak flow, 
its importance in arid climates during the dry 
season is out of the question. 

Population is increasing in the region at a 
mean annual rate of 1% for all five watersheds. 
At the present, there are currently above 10 M 
inhabitants on these watersheds and it is ex- 
pected they would reach above 12 M by the year 
2030. This new scenario envisages an increment 
on water demands of above 20% for municipal 
and potentially for industrial and irrigation 
supplies, since these sectors of the economy are 
well linked to population growth. 


Discussion 


The river flow regulation by the build up 
of reservoirs, tapping of deep groundwater 
resources, municipal discharge, and irrigation 
excess into streams is consistently noted by 
the contrasting effect of increased base flows 
and reduced direct flows of 7% of the studied 
hydrometric stations. Large reservoirs had been 
constructed having an important footprint on 
river flow. The storage capacity has increased 
exponentially in the last four decades and 
currently Northern Mexico has an ordinary 
storage capacity of above 30 000 Mm, placed 
in 10 major reservoirs that each has a maximum 
ordinary storage capacity larger than 1 000 Mm* 
in all five watersheds. 

The management of natural resources 
partially contributes to disturb river flow com- 


ponents as well. A coupled effect of increased 
direct and reduced base flow was observed in 
several gauging stations. Although there must 
be other causes that explain this contrasting 
effect, deforestation and degradation of soil re- 
sources appear to play an important role in the 
promotion of processes leading to reduced base 
flows (groundwater recharge) and increased 
direct flows (reduced infiltration and increased 
Hortonian surface runoff). Land-use changes 
from native forests to agriculture or grasslands 
lead to increased surface runoff and soil ero- 
sion rates (Návar éz Synnott, 2000a; 2000b). 
Therefore, the disappearance of Tamaulipan 
thornscrub forests that correlates well with the 
expansion of intensive agriculture and mining 
in Mexico's northeastern plains may in turn 
reduce groundwater recharge. 

Deforestation modifies the components of 
the hydrological cycle and river discharge shows 
well its fingerprint. The annual deforestation 
rate of Mexico's northeastern ecosystems ap- 
proximates 1.7% caused primarily by intensive 
agriculture and gas mining (Návar-Cháidez, 
2008). The deforestation rate approximates 
1.70% in the Pacific Ocean coastal Plains, in 
the RH 10 and 11, in agreement with land-use 
change studies conducted on tropical dry for- 
ests. Temperate forests distribute on the high 
mountain ranges and have a mean annual de- 
forestation rate of 0.50% (Ramos-Juárez, 2005). 
In these ecosystems, deforestation practices 
include overgrazing, forest wildfires, pests and 
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diseases, illegal timber harvesting, unplanned 
construction of forest roads, among others. The 
livestock carrying capacity is sometimes be- 
tween 2 to 5 times higher than the recommended 
one (Manzano € Návar, 2000; Manzano, Návar, 
Pando, € Martínez, 2000). 

Deforestation practices may also further 
magnify anomalies of river flow regimes in the 
near future, since projections under conven- 
tional scenarios of the area, previously covered 
by Tamaulipan thornscrub up to 80 490 km? 
in the 1950's (Návar-Cháidez, 2008) would be 
reduced to 22 000 km? by the year 2020 and 4 974 
km? by the year 2100. Temperate forests would 
be reduced 10% for the year 2020 and to 40% 
by the year 2100 of the total forest area mea- 
sured for the year 2000. Deforestation practices 
are expected to further reduce base flow and 
increase direct and peak flows of most upland 
watersheds. 

Climate change may also partially magnify 
river flow oscillations. However this control of 
river discharge variability cannot be properly 
quantified and we must rely on coarse projec- 
tions. Mulholland et al. (1997) and IPCC (2014) 
using global circulation models predicted 
northeastern Mexico would potentially receive 
less rainfall (between -5 to -20%) and produce 
less river flow (between 5-25% less discharge) 
because of increasing global temperatures. 
Méndez-González, Návar-Cháidez and 
González-Ontiveros (2008) have noted statistical 
evidence that the number of rainy days and total 
annual rainfall are steadily declining in the last 
four decades in several climatic stations across 
northeastern of Mexico. 

Currently, climate variability is the most 
prominent source of river flow anomalies since 
river flow parameters are well correlated with 
large scale climate events. Návar-Cháidez and 
Lizárraga-Mendiola (2013) associated large- 
scale climate events to discharge variables for a 
Mexico's northern forested watershed; e.g., the 
Southern Oscillation Index, SOI was negatively 
correlated, the PDO was positively correlated 
and the AMO was negatively correlated with 
river discharge. 


Given this state of understanding, strength- 
ening sustainable management practices of wa- 
ter resources should be addressed in parallel to 
tapping new sources of river flow in watersheds 
with high erratic river discharge, high popula- 
tion and industrial growth. Rivers, lakes and 
aquifers must supply all conventional water 
(public, household, agricultural, and industrial) 
and environmental (environmental flow for 
ecosystems). Because of the erratic frequency 
of climatic events such as hurricanes, cyclones, 
and tropical depressions, storage in constructed 
reservoirs must be increased to collect and store 
most river flow in times of high availability. In 
watersheds characterized by long term dry 
spells and important regional development 
rates, when storages are depleted, other sources 
of water supply such as the inter-basin transfer 
must be sustainably tapped. River flow and 
groundwater recharge must also be protected 
in receiving and donor watersheds (Postel, 2000; 
Schmanadt et al., 2000) and the environmental 
flow must be planned, conveyed and enforced 
in regulated rivers (Návar-Cháidez, 2014). 

Productivity must also be doubled for every 
drop of water diverted from reservoirs, streams, 
and aquifers in order to meet conventional 
supplies as well as to conserve river flow for 
meeting environmental conservation (Postel, 
2000). The current situation and likely future 
deem scenarios envisages the need for the 
implementation of the following sustainable 
practices to enhance sustainable management 
of water resources: (a) increasing storage by the 
construction of reservoirs, proper aquifer man- 
agement, and other strategies such as the imple- 
mentation of the environmental flow regime of 
regulated rivers, (b) reducing per capita water 
use, (c) increasing the efficiency of water use in 
the municipal, industrial, and the agricultural 
sector, (d) treating all raw sewages before they 
are discharged into streams and reservoirs with 
the objective to avoid the impairment of other 
sources of water supply and to be reused further 
below the watershed, (e) valuing correctly water 
resources, (f) among others. 
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Conclusions 


River flow is in transition in several of Mexico's 
northern watersheds. Currently, base and direct 
flow account for approximately 50% each of 
the total daily river discharge. A little over 40% 
of each of five river flow parameters presents 
temporal tendencies with statistical significance 
for the period of 1940 to 1999. On the average 
26% of each of the five parameters evaluated 
had a negative tendency stressing the river flow 
decline over time. Climate variability, river flow 
regulation, deforestation, land degradation, 
and potentially subtle climatic changes appear 
to explain the downward temporal monotonic 
tendencies. It is important to strengthen sus- 
tainable management practices of river flow, 
parallel to sustainable tapping new sources 
of river flow to meet conventional as well as 
environmental supplies in watersheds with high 
regional demographic and economic growth, 
with the objective of doubling productivity for 
each drop of water diverted from streams, lakes 
and aquifers. 
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